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1 .  INTRODUCTION 


The  purpose  of  the  reported  investigation  is  to  provide  insight  into 
liquid  monopropellant  jet  behavior.  More  specifically,  the  mechanisms  for 
droplet  formation  anticipated  in  the  combustion  chamber  of  existing  gun 
designs  are  analyzed. 

Emphasis  is  on  hydrodynamic  considerations;  i.e.  combustion  of  the 
liquid  is  not  Included.  While  this  is  certainly  a  simplification  of  a  very 
complex  process,  it  should  provide  insight  into  the  details  of  an  Important 
consideration  for  droplet  formation  in  a  relatively  straightforward  manner. 

The  following  chapter  provides  a  literature  survey  and  discussion  of 
possible  mechanisms  for  droplet  formation.  Then  stability  is  analyzed 
focusing  on  aerodynamic  interaction  as  the  major  mechanism  for  study.  The 
results  are  applied  to  a  typical  set  of  gun  data  and  a  comparison  is  made 
to  a  parallel  analysis  neglecting  aerodynamic  interaction  but  assuming 
a  uniform  burning  rate.  A  final  chapter  compares  the  results  with  trends 
indicated  by  other  analyses  and  actual  experimental  findings.  Extensions 
of  the  present  work  that  would  provide  further  insight  are  discussed. 
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2.  MECHANISMS  FOR  DROPLET  FORMATION  RELEVANT  TO 
THE  DESIGN  OF  LIQUID  MONOPROPELLANT  GUNS 

2.1  Introduction 

As  mentioned  earlier,  the  focus  of  this  study  will  be  on  droplet  formation 
from  a  strictly  hydrodynamic  viewpoint;  i.e.  combustion  will  not  be  included. 
Therefore,  droplet  mechanisms  that  are  applicable  to  any  liquid  jet  can  be 
discussed  and,  as  expected,  much  information  is  available,  ranging  from  simple 
hydrodynamic  theory  on  falling  drops  to  spray  development  in  diesel  engines 
or  agricultural  applications.  However,  very  little  analysis  is  available  that 
pertains  directly  to  the  design  of  liquid  monopropellant  guns.  In  the  material 
that  follows,  a  description  of  the  basic  breakup  regimes  will  be  given  and 
mechanisms  for  droplet  formation  discussed.  Then  typical  gun  data  will  be 
used  to  identify  the  regions  of  most  interest,  and  various  mechanisms  that 
are  applicable  will  be  summarized. 

2.2  Jet  Breakup  Regimes 

Jet  breakup  is  usually  discussed  in  reference  to  a  diagram  of  jet  length 
vs.  velocity.  Figure  1,  accompanied  by  a  visual  description  such  as  Figure  2. 

The  'Breakup  Length'  in  Figure  1  is  the  coherent  portion  of  the  liquid  jet,  or 
unbroken  length.  These  are  simplified  pictures  and  there  is  even  some  disagree¬ 
ment  as  to  the  shape  of  Figure  I  in  regions  of  high  velocity.  It  is  also  important 
that  Figure  1  is  qualitative  in  nature;  i.e.  it  is  usually  drawn  for  a  certain 
diameter  jet  with  given  jet  and  surrounding  properties  [typically  the  jet  is 
a  liquid  (water)  and  the  surroundings  a  gas  (air) ] . 

The  low  velocity  regions  of  the  jet  have  been  studied  in  much  detail, 
since  they  arc  amenable  to  mathematical  analysis.  In  the  region  of  very  small 
velocity,  'A'  on  Figures  I  and  2,  the  liquid  simply  drips  out  of  the  nozzle 
and  does  not  form  a  wave  at  all.  The  linear  region  'B'  (Figures  1  and  2) 

is  called  the  Rayleigh  region,  and  was  first  described  in  Rayleigh's  early 

* 

studies  on  wave  instability  [I].  Breakup  is  due  to  the  unstable  growth  of 
*  References  are  listed  on  page  23. 
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axisymmetric  surface  waves  caused  by  surface  tension,  resulting  In  droplets 
with  sizes  larger  than  the  jet  diameter.  Instability  theory  (discussed  in 
more  detail  in  the  next  chapter  with  regard  to  our  present  application) 
gives  the  wavelength  for  separation  as: 

X  -  9.02a  (1) 

where  the  wavelength  X  is  proportional  to  the  droplet  diameter  d.  (Symbols 
are  defined  in  the  Nomenclature) .  The  maximum  growth  rate  is 


u>m  -  0.118 


p!  a3 


The  time  for  droplet  separation  is  estimated  as: 


tm  -  “  -  8.47- 
in 


The  function  describing  the  linear  region  was  arrived  at  by  the  theory  of 
Weber  and  demonstrated  in  the  experiments  of  Haenlein  [2],  This  analysis 
includes  the  effects  of  viscosity  as  well  as  surface  tension. 


£.12[(^)s+3(^]; 


The  first  dimensionless  group  on  the  r.h.s.  of  the  above  equation  is  the 


Weber  number 


2pi  U02a 


which  is  a  ratio  of  the  momentum  force  4pj  UQ2  a^  to  the  surface  tension 
force  (2oa).  This  particular  parameter  can  be  based  on  either  the  liquid 
density  (as  above)  or  the  gas  density,  and  is  frequently  associated  with  both 
initial  jet  breakup  and  secondary  droplet  breakup.  Large  Weber  numbers  would, 
for  example,  indicate  that  momentum  forces  were  much  larger  than  surface 
tension  forces,  and  hence  breakup  would  likely  occur.  However,  there  is 
no  "critical"  Weber  number  that  will  successfully  predict  breakup. 

The  second  dimensionless  group  on  the  r.h.s.  of  relationship  (4)  is  the 
ratio  of  the  Weber  number  to  the  Reynolds  number. 
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The  Reynolds  number  Is  the  ratio  of  the  momentum  force  4p}U0^a2  to  the 
viscous  force  2yj  UD  a  and  Is  generally  used  to  characterize  laminar  and 
turbulent  flow  of  viscous  fluids.  The  ratio  of  the  two  is  then  characteristic 
of  the  ratio  of  viscous  forces  (2y ^  U0  a)  to  surface  tension  forces  (2oa). 

Since  viscous  forces  are  always  important  for  high  speed  jet  breakup,  this 
parameter  will  be  important  in  our  study.  High  values  will  Indicate  dominance 
of  viscous  forces  over  surface  tension  forces  and  vice  versa. 

Weber's  contributions  to  the  linear  breakup  region  also  provide  corrections 
to  relationships  (1),  (2),  and  (3)  which  account  for  the  effect  of  fluid 
viscosity.  The  wavelength  of  separation  is  given  by: 

3u  1 


\  «  8.89  a  1  + 


L  V  2a  Pi  o  J 

The  maximum  growth  rate  is  given  by: 

6yi  a  ] 
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and  the  time  for  droplet  separation  is  given  by: 

6ui  a 


(8) 
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(9) 


In  a  further  analysis,  Weber  included  the  pressure  effects  of  an  inviscid 
surrounding  [3].  This  predicted  a  maximum  in  the  breakup  region  and  the 
existence  of  region  C  (see  Figures  1  and  2) .  Weber  indicated  that  aerodynamic 
forces  tend  to  propagate  both  symmetric  and  transverse  disturbances,  leading 
to  shorter  breakup  lengths.  The  regime  in  which  these  aerodynamic  forces 
tend  to  dominate  is  called  the  first  wind-induced  regime.  This  regime  was 
also  studied  in  more  detail  by  Grant  and  Middleman  [4].  They  altered  Weber's 
theory  to  produce  a  better  fit  to  experimental  data  and  also  studied  the 
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turbulent  region,  shown  as  'D'  in  Figures  1  and  2.  In  the  turbulent  region, 
viscous  forces  tend  to  stabilize  the  jet  (in  terms  of  breakup  length).  Here 
breakup  occurs  in  the  form  of  ligaments  of  various  sizes. 

As  the  velocity  is  further  increased,  theory  becomes  less  appropriate, 
and  the  jet  definition  less  clear.  Eventually  atomization  becomes  more 
dominant  and  the  jet  breakup  length  is  reduced.  The  word  "atomization"  was 
probably  introduced  in  1875  by  Commander  Isherwood  of  American  Navy,  who  made 
a  study  of  oil  burning  for  naval  purposes  [5].  Atomization  shall  be  referred 
to  here  as  the  formation  of  small  droplets  taking  the  shape  of  a  conical 
spray  (in  the  case  of  a  round  jet).  This  formation  begins  at  the  end  of  the 
jet  in  region  D  and  approaches  the  outlet  of  the  nozzle  as  the  jet  length 
is  reduced  in  region  E,  called  the  second  wind- induced  regime  in  Figure  1. 

At  full  atomization  (region  F) ,  the  spray  begins  at  the  nozzle  exit,  although 
there  is  evidence  that  a  liquid  core  still  exists.  The  terminology  "primary" 
as  well  as  "secondary"  atomization  is  also  sometimes  used.  The  former  refers 
to  the  droplet  formation  process  (the  focus  of  this  study),  and  the  latter 
refers  to  the  disintegration  of  previously  formed  droplets. 

2.3  Additional  Factors  Contributing  to  Jet  Breakup 

The  representation  given  by  Figure  1  is  much  too  simplified.  An  alternate 
way  to  look  at  flow  regimes  is  to  consider  the  Ohnesorge  number: 


This  relationship  is  the  square  root  of  the  Weber  number  divided  by  the 
Reynolds  number  and  hence  is  Independent  of  the  velocity  but  includes  both 
surface  tension  and  viscosity.  Some  jet  regime  boundaries  are  indicated  by 
a  plot  of  this  relationship  vs.  Reynolds  number  as  in  Figure  3  [6].  The 
Ohnesorge  number  was  also  used  by  Grant  and  Middleman  to  predict  the  maximum 
in  the  laminar  breakup  curve  [A]. 


In  Figure  3,  the  solid  straight  line  to  the  Immediate  right  of  the  region 
denoted  as  the  Rayleigh  regime  Is  the  beginning  of  vhe  first  wind-induced 
regime,  region  C  on  Figure  1.  The  lines  marked  "Ohnesorge"  and  "Miesse" 
indicate  the  beginning  of  atomization.  The  word  "atomization"  is  used  in  the 
literature  to  denote  several  classes  of  jet  breakup,  and  it  is  not  clear 
whether  the  intent  here  is  to  describe  a  boundary  for  region  E  or  F  on  Figure  1. 
The  dashed  lines  are  intended  to  show  the  effect  of  gas  density.  There  is 
a  shift  in  the  wind— induced  and  atomization  boundary  to  lower  Reynolds  numbers 
as  the  gas  is  compressed.  (Other  lines  are  for  atmospheric  pressure). 

The  effect  of  gas  density  has  been  discussed  by  others  in  more  detail. 

Ranz  [7]  characterized  sprays  by  using  the  Weber  number  based  on  gas  density, 
among  other  parameters,  where: 


(ID 


Sterling  and  Sleicher  [8]  also  proposed  a  relation  between  the  Weber  number 
and  Ohnesorge  number  as  a  boundary  between  the  Rayleigh  and  first  wind- 
induced  regime.  Reitz  [9]  has  summarized  these  effects  in  Figure  4  where 
regime  boundaries  are  Indicated  as  a  function  of  density  ratio,  Weber  number, 
and  velocity.  As  expected,  atomization  is  more  likely  to  occur  as  gas 
density  increases  due  to  the  stronger  aerodynamic  forces  (momentum  forces 
of  the  gas) . 

The  curves  on  Figure  4,  as  on  Figure  3,' indicate  proposed  boundaries  between 
the  Rayleigh,  wind- induced,  and  atomization  regimes.  (Note  the  gas  Weber  number 
is  given  by  Weg  rather  than  We2) .  Ranz  [7]  proposed  that  wind-induced  jet  breakup 
(presumably  regions  C  and  D  on  Figure  1)  are  bounded  by 
.4  <  We £  <  13 
with  atomization  where 


We  2  >  13 
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The  boundary  between  the  Rayleigh  regime  and  first  wind-induced  regime  (Regions  B 
and  C  in  Figure  1)  was  proposed  by  Sterling  and  Streicher  [8]  to  be  given  by 
We 2  -  1.2  +  3.4  Oh*9 

(Note  that  Oh  is  given  by  Z  on  Figure  4).  Miesse  [6]  proposed  the  onset  of 
atomization  (beginning  of  region  E)  to  be  given  by 
We 2  -  40.3. 

Another  important  factor  influencing  jet  breakup  is  nozzle  design,  and 
the  easiest  parameter  to  quantify  in  this  regard  is  the  length-to-diameter 
ratio  (L/2a).  Reitz  and  Bracco  [10]  have  reported  results  of  this  effect, 
as  well  a§  that  of  density  ratio,  on  jet  breakup.  They  studied  several 
nozzle  designs,  operating  conditions,  and  working  fluids,  and  arrived  at 
semi-empirical  correlations  for  spray  angle,  core  length,  and  droplet  size. 

The  spray  angle  for  the  atomization  region  is  given  by: 


The  function  'f'  is  found  from  stability  theory  and  will  be  discussed 
in  the  next  chapter.  It  Is  illustrated  on  Figure  11.  The  constant  'A'  is 
determined  empirically.  For  their  range  of  experiments,  Reitz  and  Bracco 
[10]  found  A  to  be  givn  by  the  equation: 


A 


3  + 


L/2a 

3.6 


(13) 


Data  are  presented  for  spray  angle  and  two  limiting  cases  considered. 
In  the  first  case,  the  density  ratio  (liquid  to  gas)  is  high  and  the 
viscous  forces  are  small;  i.e. 
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For  this  situation.  Figure  5  shows  spray  angle  as  a  function  of  nozzle  design 
and  density  ratio.  The  design  curves  represent  equation  (12)  with  the  value 
of  A  chosen  to  provide  the  best  fit  for  all  the  data  for  a  particular  nozzle. 
Only  a  few  data  points  are  shown.  The  "goodness  of  fit"  is  a  measure  of  the 
validity  of  equation  (12).  The  'K'  curve  is  a  semi-empirical  relation  showing 
the  demarcation  between  atomization  at  the  nozzle  and  atomization  downstream 
of  the  nozzle.  In  Figure  1,  this  corresponds  to  the  division  between  region  E 
and  region  F.  The  region  above  the  'K'  curve  represents  fully  atomized  flow. 
As  indicated  before  (see  Figure  4),  as  the  gas  density  approaches  that  of  the 
liquid,  the  jet  is  more  likely  to  fully  atomize. 

The  second  limiting  case  is  for  more  viscous  flows  where 


<  <  1 


P2  \  Pi 


These  results  are  found  in  Figure  6.  Once  again,  a  'K'  curve  is  given  and 
the  region  above  the  curve  represents  full  atomization.  As  mentioned  earlier, 
viscosity  tends  to  stabilize  the  jet. 

Note  that  the  nozzle  design  curves  indicate  that  equation  (12)  is  not  as 
accurate  for  the  more  viscous  flows  (Figure  6  as  compared  to  Figure  5) . 

There  are  other  factors  that  influence  jet  behavior  in  addition  to  the 

fluid  properties,  aerodynamic  effects,  and  nozzle  geometry  given  by  the  L/2a 

ratio.  These  are  discussed  in  detail  by  Reitz  [9]  and  are  related  to  nozzle 

design,  fluid  supply  conditions,  and  the  physical  rearrangement  of  the  flow 

as  it  leaves  the  nozzle.  If  the  nozzle  has  a  sharp  inlet,  the  flow  may  not 

attach  to  the  nozzle  wall.  This  gives  a  more  stable  jet  than  one  which 

reattaches,  since  the  latter  is  more  turbulent  due  to  cavitation.  If  the 

flow  leaving  the  nozzle  is  laminar,  its  profile  is,  of  course,  less  uniform 

than  if  it  were  turbulent.  It  is  postulated  that  the  redistribution  of 

energy  as  the  flow  adjusts  to  the  removal  of  the  nozzle  produces  radial 

velocities  that  contribute  to  jet  breakup.  There  is  some  experimental 

evidence  that  fully  developed  turbulent  jets  (in  the  nozzle)  are  more 

stable  than  laminar  ones.  In  any  flow,  the  abrupt  change  in  boundary 
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conditions  leaving  the  nozzle  affects  the  tangential  stresses  and  could 
produce  short  wavelength  surface  waves  that  promote  atomization.  Finally, 
liquid  supply  pressure  oscillations,  which  are  commonly  found  In  practical 
fuel  injection  systems  can  also  have  an  effect  on  jet  breakup. 

2.4  Factors  Affecting  Droplet  Formation  In  Liquid  Monopropellant  Guns 


Many  of  the  above  mechanisms  relate  to  droplet  formation  in  liquid 
monopropellant  guns.  To  Isolate  the  problem  somewhat,  several  key  variables 
can  be  computed  for  "typical"  gun  data.  In  this  stage  of  development, 
"typical"  refers  to  test  data  which  are  currently  being  compared  to  various 
analytical  models.  Data  are  given  in  Table  1  for  the  first  few  milliseconds 
of  such  a  test.  Parameters  computed  include  L/2a,  the  Reynolds  number, 

Weber  numbers,  the  Ohnesorge  number,  the  density  ratio,  and  the  atomization 
parameter  in  equations  (14)  and  (15);  l.e. 


fl  | -a- 

P2  \ 


For  comparison  purposes,  figures  previously  presented  all  represent  particular 
fluids  and  nozzle  designs  but,  of  course,  those  with  dimensionless  groups  are 
somewhat  more  useful  than  others.  The  high  Reynolds  numbers  (>  10^)  and  low 
Ohnesorge  numbers  indicate  the  probability  of  atomization,  upon  referral  to 
Figure  3.  Even  though  Figure  4  is  for  water  jets  into  air,  the  extremely 
large  Weber  numbers  indicate  liquid  breakup  is  probable.  The  values  of  the 
density  ratios  and  abscissa  of  Figure  6  (also  tabulated)  indicate  possible 
spray  formation  when  compared  to  Figure  5.  The  L/2a  ratios  are  similar 
but  the  specific  nozzle  designs  are  quite  different.  Also,  although  spray 
formation  is  evident,  it  is  not  clear  whether  it  begins  at  the  nozzle  exit 
or  further  downstream. 

The  data  do  support  evidence  of  atomization  and  hence  legitimize  further 
study.  The  analysis  that  follows  does  consider  only  aerodynamic  effects  and 
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is  applicable  to  fluid  breakup  anywhere  this  mechanism  is  dominant.  Since 
we  are  clearly  not  in  a  regime  dominated  by  capillary  forces,  this  seems 
feasible.  However,  it  is  quite  clear  that  mechanisms  due  to  nozzle  design, 
fluid  supply  conditions,  and  the  physical  rearrangement  of  the  flow  are 
not  considered.  The  first  two  of  these  are  simply  not  amenable  to  analysis, 
since  there  is  no  experimental  or  analytical  data  on  flow  inside  the  nozzle 
itself.  The  third  could  be  considered,  but  such  rearrangement  would  augment 
droplet  development.  Preliminary  indications  are  that  the  aerodynamic 
theory  overestimates  the  amount  of  flow  in  droplet  form,  and  inclusion  of 
this  mechanism  would  only  exacerbate  this  situation;  i.e.  any  additional 
mechanism  considered  beyond  that  in  the  next  chapter  should  be  flow 
stabilizing  and  not  destabilizing. 


I 

i 
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3 .  ANALYSIS 

In  Che  analysis  Chat  follows,  it  will  be  assumed  Chat  breakup  occurs 
because  of  instabilities  developed  in  a  liquid  jet  caused  by  aerodynamic 
forces.  The  procedure  is  outlined  by  Reitz  [9  ]  and  Reitz  &  Bracco  [10] 
with  details  given  in  the  analyses  of  Levich  [11]  and  Taylor  [12].  For 
the  analysis  an  axisymmetric  cylindrical  jet  will  be  considered.  Results 
will  be  applied  to  a  geometry  resembling  that  of  a  liquid  monopropellant  gun. 

The  linearized  Navier-Stokes  equations  for  the  small  axisymmetric 
fluctuating  (perturbation)  velocities  u^,  vj,  and  pressures  pj  (i  ■  1  for 
liquid,  i  -  2  for  gas)  are: 


ir  +  £  £  (rvi>  ■ 0> 

But 

+  Ui(r)  +  v 
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Symbols  are  further  defined  in  the  Nomenclature.  Assuming  the  surface  wave 
elevation  n  <  <  a,  the  boundary  conditions  to  a  first  order  approximation 
can  be  written  as  (see  Figure  7): 
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Solutions  are  obtained  by  Introducing  a  velocity  potential  $,  and 
stream  functions  4>i»  ^2  the  form 


♦x  ■  + 

(7) 

*1  -  »1(r).lkz  +  “t. 

(8) 

*2  “  [U2 Cr)  -  l(u/k)]  C3  .ikz  +  wt  f(r) 

(9) 

Substitution  of  (7)  and  (8)  Into  equations  1-3  for  the  liquid  component 
yields  a  solution  in  the  form  of  modified  Bessel  functions  In: 

*1  "  C1I€[(kr),  (10) 

-  CjrliClr),  (11) 

l2  ■  k2  +  u/vj  *  (12) 


For  the  gas  component,  an  Orr-Summer field  equation  for  f(r)  Is  solved  to 

give  the  gas  pressure: 

/t»  u  ^o  (k8)  ikz  +  «t 

p2  “  "  p2(uo  -  k  Kj(ka)  n°e  » 

where  the  Kq  are  Bessel  functions  of  the  second  kind. 


(13) 


The  constants  Cj_,  C2»  C3  are  solved  for  by  considering  the  boundary 
conditions  equations  4-6.  This  results  in  the  following  dispersion 


relationship: 
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where  the  prime  denotes  differentiation.  The  above  equation  relates  growth 
rate  «  to  wavelength  X,  but  is  not  easily  solved.  A  limiting  solution  for 
ka  •  (small  wavelength  disturbances)  where  the  droplet  size,  r  ^  X,  is  small 
is  applicable  to  atomization.  In  this  case  the  Bessel  functions  in  (14)  can 
be  replaced  by  their  asymptotic  values,  and  (14)  becomes: 
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(u  +  2vjk2)^  +  ak^/pj  -  4v^  k^(fc2  +  oj/vj)^ 


(15) 


+  (u  +  iU0k)2(p2/p1)  -  0 

This  equation  was  arrived  at  independently  by  Levich  (see  his  equation  125.17) 
for  the  case  of  high  velocity  jet  breakup  and  short  wavelengths,  and  by 
Taylor  (see  his  equation  16)  for  ripples  formed  on  an  Infinite  viscous  liquid. 
This  is  expected,  since  for  kr  -*■  •  curvature  effects  would  be  unimportant; 
i.e.  note  (15)  above  is  Independent  of  r.  Taylor  considered  a  further 
simplification  assuming: 


P 2/Pi  <  <  1  (16) 

In  our  case,  P2/P1  ^  .014  to  .13.  Taylor  separated  equation  (15)  above 
into  real  and  imaginary  components.  With  assumption  (16),  he  obtained  for 


the  real  part  of  the  disturbance  growth  rate,  to: 
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where  B  is  the  droplet  formation  parameter: 
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x  is  the  wavelength  parameter: 

P2Uo2 
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(i:) 


and  g  is  given  as  a  function  of  x  on  Figures  8  and  9. 

Note  that  the  droplet  formation  parameter  has  also  been  used  to  predict 
the  onset  of  atomization  (see  equations  14  and  15  in  Chapter  2). 

Waves  identified  with  a  growth  rate  described  by  relation  (17)  are 
unstable.  Separation  into  drops  occurs  for  a  maximum  value  of  g;  i.e. 

g(B,x)  -  gm(B,xm)  (20) 


where  gm  and  Xq  define  the  values  of  u  and  k  for  the  wave  that  forms  a  drop. 
Figures  8,  9,10,  and  11  show  the  location  of  these  maxima.  Note  the  shift 
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to  longer  wavelengths  (larger  x)  with  smaller  values  of  B.  Increasing  x  means 
decreasing  k  and  Increasing  wavelength  since: 

k  -  2ir/X  (21) 

This  shows  that  highly  viscous  flow,  giving  smaller  B,  leads  to  longer  waves 
and  less  viscous  flow  produces  shorter  waves.  In  our  case  B  can  range  from 
an  order  of  magnitude  of  20  to  .001  with  much  of  the  data  having  values  less 
than  1 . 

A  wave  of  small  amplitude  grows  on  the  surface  until  the  crest  is  detached 
and  a  drop  is  formed.  The  diameter  of  the  drop  can  be  associated  with  the 
wavelength  of  this  detached  wave: 


d  ■  C  1, 


In  this  analysis,  the  simplification  will  be  made  that  the  constant  C  *  1. 
As  mentioned  above,  Xm  is  associated  with  x„,  since: 

Xm  *  2  (. 

The  growth  rate  ui  is  used  to  estimate  the  amount  of  time  required  to  form 
a  drop.  During  a  time: 


the  maximum  wave  grows  by  an  amount 

.  e 

A  reasonable  estimate  for  the  time  to  form  a  drop  is  a  value  proportional 
to  tn  and,  in  this  analysis,  the  proportionality  constant  will  be  considered 


unity. 


The  mass  of  fluid  detached  from  unit  area  of  the  surface  can  be 


written  as 


Pj  it  D  d  l  Xffl 


P,  Xm 
1  ® 


ir  D  d  l 


(25) 


where  di  is  a  differential  length  of  the  jet  and  D  is  the  diameter  of  the 
jet  (in  our  case  the  diameter  of  the  bolt  plus  twice  the  jet  thickness). 
The  mass  detached  per  unit  area  per  unit  time  is: 


Pl^ 


(26) 


Combining  equations  (19)  and  (21)  yields: 

2r  2  it  a*,, 

^  ^  "  P2  U02 

where  the  subscripts  m  again  refer  to  maximum  values;  i.e.  the  values 
corresponding  to  drop  formation.  Equations  (17)  and  (24)  combine  to  get 
the  time  for  drop  formation. 


(27) 


2  u03  r  p23  l !  1  -1 
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Equations  (26) ,  (27) ,  and  (28)  combine  to  give  the  detached  mass  per  unit 
time,  per  unit  area,  as: 


m  -  4xpi 


where 


8m  *m) 


(29) 


(30) 


Figure  11  shows  fm  as  a  function  of  B. 
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4.  APPLICATION  OF  THE  ANALYSIS  TO  LIQUID  MONOPROPELLANT  JETS 


The  geometry  of  the  regenerative  liquid  propellant  gun  is  shown  in 
Figure  12.  Upon  firing,  the  annular  piston  compresses  the  fluid  in  the 
propellant  chamber  which  is  then  released  into  the  combustion  chamber  via 
the  annular  gap  created  between  the  piston  and  the  stationary  control  rod, 
or  bolt.  The  jet  is  annular  in  shape  and  is  partially  adjoined  on  its 
Inside  surface  by  the  bolt. 

4.1  Reduction  of  the  Analysis  to  Computer  Code 


From  the  previous  chapter,  the  relationship  for  mass  removed  from  the 
jet  per  unit  surface  area  and  time  is:  [see  equation  (29)  from  Chapter  3)]. 


m  •  4  r 


51  u°  f  si  )  ‘  £» 


(1) 


At  any  instant  in  time,  there  is  a  certain  surface  area  of  liquid  remaining 
intact  in  the  combustion  chamber. 

In  actual  determination  of  the  mass  removed,  two  models  were  used. 

The  first  was  discrete  in  formulation.  Over  an  increment  of  time,  dt,  the 
jet  length  is: 


di.  -  U0  dt 


and  the  surface  area  is: 


2  it  DdJL 

To  get  the  mass  removed  during  time  dt,  combining  equations  (1)  through 
(3)  yields: 

1 


(2) 


(3) 


[P2  *  1 

■  "  [4ffpl  Uo  ^  fm  J 


2 w  DdZ  dt 


(4) 


During  the  next  time  increment,  this  amount  of  mass  is  again  reduced  (the 
velocity  used  in  equation  (4)  is  now  the  velocity  at  the  present 
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time).  Also,  a  new  amount  of  mass  has  entered  and  part  of  the  new  amount 
has  been  removed.  The  mass  removed  is  distributed  into  droplets  with 

diameter  given  by  relation  (27)  in  Chapter  3,  or: 

,  2*  ax* 

”  “  PTV"  (5) 

where  the  values  of  Xq  and  UQ  are  evaluated  at  the  current  time.  An  inventory 
gives  the  mass  removed,  the  number  and  size  of  droplets,  and  the  liquid 
remaining  in  the  Intact  core  at  any  time. 

The  second  model  used  to  determine  the  mass  removed  was  continuous  in 
nature.  Equation  (1)  above  was  multiplied  by  an  appropriate  surface  area 
and  Integrated  using  an  o.d.e.  solver.  The  appropriate  area  was  considered 
to  be  that  of  a  hollow  frustrum  with  inner  diameter  equivalent  to  that  of 
the  bolt  and  outer  diameter  equal  to  the  bolt  plus  two  annuli.  The  height 
of  the  frustrum  then  corresponds  to  the  length  of  the  Intact  core.  Equation 
(5)  was  again  used  to  determine  the  droplet  diameter,  and  the  mass  and 
surface  area  of  droplets  removed  were  computed. 

The  computer  programs  for  both  methods  of  analysis  are  given  in  the 
Appendix.  The  discrete  method  has  more  physical  appeal,  but  the  continuous 
model  is  one  which  can  be  easily  incorporated  as  an  integral  part  of  existing 
gun  codes,  and  hence  satisfies  the  requirement  for  this  project. 

4.2  Consideration  of  Mass  Removal  by  Burning  (Combustion) 

Although  the  thrust  of  this  analysis  is  hydrodynamic  in  nature,  the 
models  mentioned  above  can  be  altered  to  consider  combustion  as  the  major 
mechanism  for  mass  removal  from  the  liquid  core.  Such  an  analysis  could  be 
used  to  test  the  magnitude  of  the  two  effects;  i.e.  if  the  mass  is  removed 
more  quickly  by  combustion,  then  perhaps  hydrodynamics  are  not  important, 
and  vice  versa.  The  discussion  also  pertains  only  to  initial  mass  removal; 
i.e.  mass  could  be  removed  hydrodynamically ,  form  drops,  and  the  drops 
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subsequently  burned.  The  combustion  model  Incorporates  a  simple  burning 
rate  relation  [13]: 

.  103 

Burning  Rate  (cm/s)  ■  1.64  (Pressure  In  MPa)*  (6) 

The  above  can  be  applied  to  the  surface  area  of  the  liquid  jet,  and  the 
amount  of  mass  removed  by  burning  can  then  be  determined  at  any  time. 

4.3  Application  of  the  Analysis  to  Actual  Experimental  Gun  Conditions 

The  analysis  has  been  applied  to  data  computed  from  an  existing  lumped 
model  code  which  predicts  relative  liquid  velocity,  area  of  the  annulus, 

liquid  and  gas  densities,  and  pressure  as  a  function  of  time  [13].  This 

information  provides  an  input  to  the  computer  codes  given  in  the  Appendix. 
Tables  2  and  3  illustrate  the  results  for  approximately  the  first  2.5 
milliseconds  of  liquid  entry  into  the  combustion  chamber.  In  the  table, 
the  time  is  in  milliseconds,  the  mass  in  grams,  the  length  in  cm,  and  the 
area  in  cm^.  Table  2,  which  illustrates  the  results  using  the  continuous 
approach,  shows  that  about  97Z  of  the  mass  is  atomized  after  about  2.5  milli¬ 
seconds,  a  short  liquid  core  of  about  .5  cm  remains,  and  about  4.6  x  10^  cm^ 

of  surface  area  is  formed  by  the  droplets.  The  discrete  approach,  shown 

in  the  first  four  columns  of  Table  3,  shows  very  similar  results.  Differences 
occur  because  of  the  nature  of  the  models.  Results  of  the  combustion 
analysis  are  also  given.  Note  that  very  little  mass  (less  than  1Z)  is 
removed  by  burning.  This  seems  to  indicate  that  mass  removal  by  aerodynamic 


forces  is  a  much  more  dominant  mechanism  than  mass  removal  by  burning. 


5.  CONCLUSIONS  OF  THE  STUDY  AND  RECOMMENDATIONS  FOR  FUTURE  ANALYSIS 


The  results  of  the  analysis,  when  applied  to  the  liquid  monopropellant 
jet,  indicate  that  the  flow  regime  is  in  the  atomization  region  and,  if 
data  from  a  gun  code  representing  relative  velocity,  annulus  area,  and 
fluid  densities  are  used  in  the  stability  analysis,  liquid  is  stripped 
from  the  jet  by  aerodynamic  forces,  presumably  forming  drops.  This 
mechanism  has  been  shown  to  be  much  more  dominant  than  the  combustion 
process  based  on  a  simple  burning  rate  model. 

Studies  at  the  Ballistics  Research  Laboratory  have  used  experimental 
data  in  an  "inverse  code".  By  neglecting  energy  losses  and  using  actual 
measurements,  the  amount  of  liquid  reacting  in  the  combustion  chamber  can 
be  evaluated.  When  this  is  compared  with  the  amount  of  liquid  that  has 
entered  the  chamber  (also  found  by  measurement  and  conservation  of  mass) , 
it  is  inferred  that  there  is  substantial  liquid  accumulation;  i.e.  unburned 
liquid  is  present  in  the  combustion  chamber.  From  this  result  one  would 
suspect  that  droplets  have  not  formed,  since  it  they  had  formed,  they  would 
presumably  have  burned  and  liquid  would  not  have  accumulated. 

This  presents  a  dilemma  which  bears  further  investigation.  Some  possible 
explanations  are  discussed  below  (there  are  no  doubt  more  that  perhaps  are 
not  as  obvious) . 

The  analysis  presented  in  Chapter  3  could  be  incorrect  when  applied  to 
the  present  situation.  There  are  several  possible  reasons  for  this.  The 
model  strictly  holds  for  round  jets,  not  jets  adjoined  by  a  wall  as  is  the 
case  here;  i.e.  the  geometry  is  not  consistent  with  the  model.  The  model 
also  is  inexact  in  other  ways,  most  noticeably  in  the  fact  that  it  is  for 
density  ratios  that  are  much  less  than  one  (a  term  containing  P2/P1  is 
actually  dropped  in  the  analysis).  Combustion  is  also  neglected.  Although 
this  was  considered  separately  in  a  simplistic  way,  no  attempt  was  made  to 
combine  the  two  effects. 


It  is  also  possible  that  the  jet  does  break  up,  but  does  not  form 
droplets  and  burn.  Ligaments  may  be  torn  from  the  jet  and  they  may  fall 
or  strike  other  portions  of  the  chamber,  perhaps  even  coalescing.  The 
analysis  also  is  not  very  exact  when  it  comes  to  droplet  size.  Some 
ligaments  could  be  longer  than  others,  and  constants  taken  in  the  theory 
to  be  unity  could,  in  fact,  have  almost  any  value.  If  there  is  not 
substantial  increase  in  surface  area,  burning  could  proceed  slowly  (the 
burning  analysis  actually  showed  this  in  the  limiting  case  of  no  droplet 
formation).  The  analysis  did  not  consider  what  happened  after  mass  was 
removed  from  the  jet;  i.e.  it  is  not  clear  what  aerodynamic  forces 
(relative  velocities)  are  present  once  material  has  left  the  surface. 

Another  possibility  for  the  discrepancy  is  the  existence  of  other  jet 
breakup  mechanisms,  such  as  those  mentioned  earlier  and  discussed  in  more 
detail  by  Reitz  [9].  If  the  jet  does  not  adhere  to  the  nozzle  walls,  it 
could  perhaps  be  more  stable  than  the  theory  predicts.  More  detailed 
analysis  of  flow  in  the  nozzle  (currently  under  investigation  by  others) 
might  clear  up  this  concern.  Consideration  of  flow  being  rearranged  as  it 
enters  the  combustion  chamber  would  perhaps  lead  to  more  destabilizing 
results,  although  the  fact  that  the  fluid  is  adjoined  by  a  surface  on  one 
side  could  suggest  further  stabilization. 

There  is  also  the  likelihood  that  the  input  data,  particularly  relative 
velocities  and  annulus  area  (positon  of  the  piston)  are  in  error.  A  complet 
set  of  experimental  data  does  not  exist  and  is  certainly  not  repeatable 
at  present.  Any  output  from  a  gun  code  will  be  questionable  until  this 
occurs  and  verification  takes  place.  "Inverse  codes"  that  use  experimental 
data  are  somewhat  more  reliable,  but  a  careful  error  analysis  of  all  data 
and  subsequent  exhaustive  computer  runs  covering  all  possibilities  are 
required  to  obtain  complete  confidence  in  results. 


A  related  explanation  is  that  unaccounted  for  energy  losses  could 
reduce  the  liquid  accumulation  predicted  by  the  inverse  code.  It  is  assumed 
that  the  pressure  rise  results  from  combustion  of  the  liquid  with  no  energy 
losses.  More  liquid  could  burn  with  energy  being  lost  due  to  heat  transfer. 

When  evaluating  the  above ,  some  possibilities  appear  more  likely  than 
others,  although  the  magnitude  of  error  cannot  be  determined  without  further 
study.  The  most  useful  exercise  would  be  to  conduct  conclusive  repetitive 
experiments  with  reliable  data.  However,  this  is  the  most  difficult  and 
expensive  suggestion.  Based  on  a  rather  exhaustive  body  of  literature, 
and  the  likelihood  that  at  least  velocities  are  well  predicted  by  the  code, 
it  seems  probable  that  aerodynamic  forces  do  strip  liquid  from  the  jet 
as  predicted.  It  is  not  at  all  clear,  however,  what  happens  next.  Perhaps 
droplets  are  not  formed  and,  if  they  are,  their  size  is  hard  to  determine 
with  any  accuracy.  Coalescence  due  to  interaction  among  particles  or 
adherence  to  surfaces  is  also  possible.  Another  possibility  which  is 
easier  to  check  is  the  existence  of  substantial  energy  losses  that  would 
reduce  the  amount  of  liquid  accumulated  in  the  inverse  code. 

There  is  clearly  much  work  to  be  done,  but  perhaps  the  above  analysis 
has  helped  to  somehow  enhance  our  understanding  of  the  breakup  of  liquid 
monopropellant  jets. 


TABLE  1.  VARIABLE  TABULATION  FOR  GUN  DATA 
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TABLE  2.  ATOMIZATION  PREDICTIONS  USING 
THE  CONTINUOUS  MODEL 


Mass 

Jet 

Mas  8 

Surface 

Time 

Left 

Length 

Atom. 

Area 

6.15 

0.0490 

0.1151 

0.0400 

18.98 

6.20 

0.0799 

0.1808 

0.0982 

44.93 

6.25 

0.0992 

0.2164 

0.1684 

74.84 

6.30 

0.1127 

0.2376 

0.2457 

106.67 

6.35 

0.1230 

0.2545 

0.3280 

139.96 

6.40 

0.1292 

0.2638 

0.4141 

174.60 

6.45 

0.1356 

0.2740 

0.5037 

210.82 

6.50 

0.1390 

0.2807 

0.5967 

248.92 

6.55 

0.1424 

0.2876 

0.6986 

298.06 

6.60 

0.1478 

0.2984 

0.8220 

382.93 

6.65 

0.1456 

0.2902 

0.9732 

540.93 

6.70 

0.1411 

0.2812 

1.1531 

830.56 

6.75 

0.1430 

0.2810 

1.3477 

1214.35 

6.80 

0.1535 

0.2903 

1.5601 

1709.12 

6.85 

0.1612 

0.2984 

1.7939 

2337.59 

6.90 

0.1747 

0.3128 

2.0506 

3115.51 

6.95 

0.1875 

0.3257 

2.3298 

4047.82 

7.00 

0.2007 

0.3359 

2.6303 

5134.89 

7.05 

0.2123 

0.3447 

2.9476 

6324.28 

7.10 

0.2257 

0.3559 

3.2701 

7485.80 

7.15 

0.2397 

0.3628 

3.5949 

8603.88 

7.20 

0.2393 

0.3491 

3.9325 

9833.57 

7.25 

0.2384 

0.3366 

4.2799 

11261.07 

7.30 

0.2456 

0.3383 

4.6200 

12615.46 

7.35 

0.2574 

0.3421 

4.9594 

13922.19 

7.40 

0.2677 

0.3484 

5.3001 

15193.36 

7.45 

0.2796 

0.3529 

5.6429 

16435.98 

7.50 

0.2924 

0.3606 

5.9885 

17652.54 

7.55 

0.3089 

0.3694 

6.3483 

18970.00 

7.60 

0.3359 

0.3878 

6.7480 

20667.38 

7.65 

0.3452 

0.3900 

7.1985 

22932.20 

7.70 

0.3593 

0.3923 

7.7030 

26055.22 

7.75 

0.3776 

0.4008 

8.2409 

29712.64 

7.80 

0.4071 

0.4160 

8.8182 

33955.67 

7.85 

0.4207 

0.4159 

9.4480 

39136.31 

7.90 

0.4455 

0.4191 

10.1270 

45496.38 

7.95 

0.4827 

0.4342 

10.8407 

52494.75 

8.00 

0.5244 

0.4519 

11.5963 

60203.61 

8.05 

0.5437 

0.4539 

12.4061 

69026.13 

8.10 

0.5843 

0.4664 

13.3149 

80860.73 

8.15 

0.6368 

0.4859 

14.3131 

95341.89 

8.20 

0.7162 

0.5138 

15.4083 

112671.40 

8.25 

0.7615 

0.5233 

16.6209 

134052.50 

8.30 

0.8366 

0.5388 

17.9474 

160528.90 

8.35 

0.8823 

0.5510 

19.3589 

190107.10 

8.40 

0.9486 

0.5900 

20.8876 

224487.50 

8.45 

0.9507 

0.5913 

22.5692 

267729.40 

8.50 

0.9443 

0.5873 

24.3774 

321690.50 

8.55 

0.9575 

0.5955 

26.2804 

384149.70 

8.60 

0.9437 

0.5869 

28.3386 

461050.10 

TABLE  3.  ATOMIZATION  AND  BURNING  RATE  PREDICTIONS 
USING  THE  DISCRETE  MODEL 


Atomization  Burning 


Mass 

Mass 

Jet 

Mass 

Mass 

Mass 

Time 

Atom. 

Left 

Length 

Burned 

Left 

In 

6.15 

0.0338 

0.0539 

0.1105 

0.0000935 

0.0876 

0.0877 

6.20 

0.0985 

0.0779 

0.2175 

0.0002776 

0.1761 

0.1764 

6.25 

0.1810 

0.0850 

0.2114 

0.0005505 

0.2655 

0.2660 

6.30 

0.2662 

0.0910 

0.2072 

0.0009105 

0.3563 

0.3572 

6.35 

0.3538 

0.0960 

0.2048 

0.0013570 

0.4485 

0.4498 

6.40 

0.4410 

0.1028 

0.3068 

0.0018899 

0.5420 

0.5439 

6.45 

0.5319 

0.1076 

0.3065 

0.0025123 

0.6370 

0.6395 

6.50 

0.6266 

0.1100 

0.3083 

0.0032268 

0.7334 

0.7366 

6.55 

0.7372 

0.1090 

0.3234 

0.0040439 

0.8422 

0.8462 

6.60 

0.8721 

0.1067 

0.2587 

0.0049862 

0.9738 

0.9788 

6.65 

1.0386 

0.0962 

0.3071 

0.0060781 

1.1288 

1.1348 

6.70 

1.2409 

0.0726 

0.1880 

0.0073420 

1.3062 

1.3135 

6.75 

1.4422 

0.0715 

0.2090 

0.0087941 

1.5049 

1.5137 

6.80 

1.6649 

0.0725 

0.2285 

0.0104519 

1.7270 

1.7374 

6.85 

1.9103 

0.0751 

0.2461 

0.0123342 

1.9731 

1.9855 

6.90 

2.1789 

0.0780 

0.2618 

0.0144587 

2.2424 

2.2569 

6.95 

2.4690 

0.0815 

0.2756 

0.0168338 

2.5337 

2.5505 

7.00 

2.7794 

0.0877 

0.2874 

0.0194695 

2.8476 

2.8671 

7.05 

3.0981 

0.0997 

0.2889 

0.0223722 

3.1754 

3.1978 

7.10 

3.4180 

0.  j.120 

0.2813 

0.0255337 

3.5044 

3.5300 

7.15 

3.7404 

0.1234 

0.2743 

0.0289488 

3.8348 

3.8633 

7.20 

4.0783 

0.1230 

0.2679 

0.0326128 

4.1687 

4.2013 

7.25 

4.4186 

0.1240 

0.2615 

0.0365214 

4.5061 

4.5426 

7.30 

4.7514 

0.1343 

0.2562 

0.0406747 

4.8450 

4.8857 

7.35 

5.0869 

0.1439 

0.2510 

0.0450702 

5.1858 

5.2308 

7.40 

5.4253 

0.1531 

0.2462 

0.0497022 

5.5286 

5.5783 

7.45 

5.7667 

0.1618 

0.2420 

0.0545678 

5.8740 

5.9285 

7.50 

6.1114 

0.1702 

0.2380 

0.0596757 

6.2219 

6.2816 

7.55 

6.4789 

0.1799 

0.2482 

0.0650272 

6.5937 

6.6588 

7.60 

6.8905 

0.1924 

0.2714 

0.0706453 

7.0122 

7.0829 

7.65 

7.3609 

0.1919 

0.2927 

0.0765591 

7.4762 

7.5528 

7.70 

7.8777 

0.1894 

0.3119 

0.0828009 

7.9842 

8.0670 

7.75 

8.4247 

0.2016 

0.3291 

0.0893850 

8.5369 

8.6263 

7.80 

9.0140 

0.2172 

0.3443 

0.0963163 

9.1348 

9.2312 

7.85 

9.6618 

0.2209 

0.3576 

0.1036224 

9.7790 

9.8826 

7.90 

10.3539 

0.2277 

0.3691 

0.1113253 

10.4703 

10.5816 

7.95 

11.0755 

0.2511 

0.3790 

0.1194201 

11.2072 

11.3266 

8.00 

11.8398 

0.2759 

0.3873 

0.1279144 

11.9878 

12.1158 

8.05 

12.6984 

0.2848 

0.4084 

0.1368446 

12.8464 

12.9832 

8. 10 

13.6709 

0.2903 

0.4408 

0.1462460 

13.8149 

13.9612 

8.15 

14.7288 

0.3176 

0.4692 

0.1561203 

14.8903 

15.0464 

8.20 

15.8860 

0.3560 

0.4938 

0.1665039 

16.0756 

16.2421 

8.25 

17.1688 

0.3791 

0.5146 

0.1774464 

17.3705 

17.5480 

8.30 

18.5547 

0.4053 

0.5321 

0.1889534 

18.7711 

18.9600 

8.35 

20.0157 

0.4621 

0.5466 

0.2010231 

20.2769 

20.4779 

8.40 

21.6219 

0.4809 

0.5717 

0.2137051 

21.8891 

22.1028 

8.45 

23.4150 

0.4192 

0.6092 

0.2270412 

23.6072 

23.8342 

8.50 

25.3274 

0.3454 

0.6469 

0.2410622 

25.4317 

25.6728 

8.55 

27.3489 

0.2860 

0.6904 

0.2558193 

27.3791 

27.6349 

8.60 

29.5691 

0.1680 

0.7397 

0.2713833 

29.4657 

29.7371 

QUALITATIVE  PICTURE  OF  JET  BREAKUP 


C 


<  •  • 


k  os 

<  3  < 
5  ^ 


=i  < 

Q.  Ui 

^  GO 

u 


i 

< 


FIRST  AERODYNAMIC 

AERODYNAMIC  REGIME 

REGIME  SECONDARY 

LAMINAR  ATOMIZATION 

TURBULENT 


LIQUID- WATER 
NOZZLE  do.0.025cm 


U0  (cm/s) 

FIGURE  4. 

EFFECT  OF  CAS  DENSITY  ON  JET  BREAKUP  REGIME  BOUNDARIES  [9] 
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FIGURE  7.  FLOW  FIELD  FOR  STABILITY  ANALYSIS 
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FIGURE  9.  DISPERSION  RELATION  FOR  SECOND  WIND  INDUCED  BREAKUP  [9 


FIGURE  10.  LOCUS  OF  MAXIMUM  POINTS  OF  FIGURES  8 


BEFORE  FIRING 


Figure  12.  A  Regenerative  Liquid  Propellant 
with  an  Annular  Piston  [13] 
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NOMENCLATURE 


C,Ci»C2*C3 


jet  radius 

constant,  equations  (12)  and  (13),  Chapter  2 

droplet  formation  parameter,  P i /p 2 C — ^rr-)^ 

U 1  Un 


constants 


drop  diameter 

annulus  outer  diameter 

function,  equation  (30),  Chapter  3 

function,  equation  (17),  Chapter  3 

disturbance  wave  number,  2tt/A 

wave  parameter,  equation  (12),  Chapter  3 

length  of  nozzle 

Intact  length  of  jet 


mass  removed  per  unit  surface  area  and  time 


pressure 


direction  perpendicular  to  axis 


fluctuating  (perturbation)  velocity  in  z  direction  (parallel  to  jet  axis) 
relative  velocity  of  the  liquid,  Uj  -  U2 
axial  velocity  profile  in  the  liquid 

axial  velocity  profile  in  the  gas 

fluctuating  (perturbation)  velocity  in  r  direction  (perpendicular  to 
jet  axis) 

velocity  of  gas  or  liquid  in  axial  z  direction 
wavelength  parameter,  p£  Uq  /ak 
direction  parallel  to  axis 


•  s’*  ,  1  W  -T 


Nomenclature  (continued) 


6  •  annulus  thickness 

X  -  disturbance  wavelength 
n  11  surface  wave  elevation  (r  direction) 
v  -  kinematic  viscosity 

V  -  absolute  viscosity 

o  ■  surface  tension 

u  -  disturbance  growth  rate  (1/time) 

4  ■  velocity  potential 

if  »  stream  function 
8  -  spray  angle 

Subscripts 

i  ■  component 

1  »  liquid  component 

2  ■  gas  component 

o  ■  jet  (at  nozzle  exit) 

m  -  maximum  value,  value  corresponding  to  drop  formation 
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APPENDIX  A 


Computer  Programs 
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o  O  O  O  o  O  •  O  C<  Ci  O  O  O  O 


Print  file  "JETC.FTN 


Pag# 


PROGRAM  TO  CALCULATE  THE  AMOUNT  OF  LIQUID  THAT  IS  ATOMIZED  IN  AN  ANNULUS 
PROGRAM  NAME  IS  JETC.FTN.THIS  VERSION  ALSO  CALCULATES  DROP  SIZE  DISTRIBUTION. 
WRITTEN  BY  NELSON  MACKEN  OCT. 86. REVISED  21  DEC  TO  ALSO  COMPUTE  BURNING  RATE. 

IN  THIS  PROGRAM, DIMENSIONS  ARE  VAR ( J ,  I ) , WHERE  J  TRACKS  LUMP  OF  MASS  THAT 
LEFT  INLET  AT  TIME  T(J),AND  I  REPRESENTS  CURRENT  TIME. 

DIMENSION  T(60) ,V(60) ,VI(60) ,VP(60) ,AV(60> ,DEL(60) , 

1  RHOL  <  60 ) , RHOG ( 60 ) ,SDR(60,60) ,SDRT<60) ,BRATE<60) , 

2  BAREA ( 60 ) , PRESS ( 60 ) , tTMBUR (0:60), TTMLEB (0:60) 

REAL  L(60) ,LNET(60) ,MIN(60) ,MNET(0:60> ,MREMO( 60,60) , 

1  MREMA (0:60,0:60), MATOM ( 60 ) , MINNET , MNETT , LNETT , MDR , 

2  LM ( 60 , 60 ) , MLEFTB (0:60,0:60) ,MBURN (60,60) 

REAL  NDRJT, NDR (60,60), NDRT ( 60 ) , NDRNET , MINTT ( 60 ) 

INTEGER  LL(60),LBURN(60> 

OPEN( 10, FILE3 'CONST. JETC' , STATUS3 1  OLD ' ) 

OPEN (20, FILE3* OUTPUT. J ETC .STATUS3 ' NEW' ) 

PI=3. 1416 
DT= .05E-3 
ST=71 .6 
VISL= .03935 
DB=1 .7 

N*51 

100  FORMAT! 7(F7.2,2X) ) 

READ( 10, 100) (T( I ) , VI ( I ) ,VP< I ) ,AV( I ) ,RHOL( I ) , RHOG ( I ) , 

1  PRESS! I), 1=1, N) 

WRITE (20, 100) (T(I),VI(I),VP(I),AV(I) ,RHQL( I ) , 

1  RHOG(I) ,PRESS(I) ,I=1,N) 

RESET  DATA  SO  USE  AVERAGES  FOR  EACH  TIMESTEP. 

N=50 

DO  170  1  =  1,  N 
T( I ) =T( 1+1 ) 

V(I)3(VI ( I)+VI ( 1+1 ) >/2. 

AV( I ) 3 ( AV( I )+AV( 1+1 ) ) /2. 

RHOL ( I ) 3 ( RHOL ( I ) +RHOL ( I + 1 ) ) / 2 . 

RHOG(I)3(RHOG(I)+RHOGC  *l))/2. 

1 70  PRESS ( I ) 3 ( PRESS ( I ) +PRE1  3 ( I + 1 ) ) / 2 . 


DO  110  1=1, N 
L( I ) *V( I ) *DT 

MIN( I ) =AV( I )*V< I )*RHOL< I )*DT 
110  CONTINUE 


C  120  FORMAT (2X,F4. 2, 2X,F6. 1 ,2X,F5.3,2X,F8.5) 

C  WRITE (20, 120) (T( I ) ,V( I ) ,DEL( I > , B ( I > ,I=1,N) 

C  i 
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Print  fit*  “JETC.FTN 


P19* 


LLL«2 

IF(LLL.EQ.l)  SO  TO  500 
DO  150  1=1, N 
150  LL ( 1 ) =0 
MINNET*0. 

SDRNET  *0 . 

NDRNET  *0 

********  PRINT  TITLE  FOR  OUTPUT  »♦**♦♦** 

PRINT*," 

PRINT*,'****  DROP  SIZE  DISTRIBUTION  ****** ' 

DO  130  1=1, N 
PRINT*,'  TIME3  ' ,T(I ) 

MREMA( I , I- 1 ) =MIN( I ) 

MNETT =0 . 

LNETT=0. 

NDRTT=0 

SDRTT=0. 

DO  140  J=1,I 
IF(LL(J).NE.l)  THEN 

B=(RH0L(J)/RH0G(J))#(ST/(VISL*V(I)))**2. 

XP=L0G10(B) 

IF1XP.LE.-5.)  F=.01*XP+.09 

IF(XP.6T.-5. .AND.XP.LE.-3. )  F«.02*XP+.14 
IF(XP.GT.-3..AND.XP.LE.-1.)  F  = . 0375*XP+ .1925 
IF(XP .GT.-l . .AND.XP.LE. 1 . )  F=.045*XP+.2 

IF(XP.GT.l)  F=.0l75*XP+.2275 

I F  <  XP . LE . -2 . )  XM=10.**<-.29*XP+.22> 

IF(XP.GT. -2. . AND.XP.LE. -1 . )  XM=10.**( - .24*XP*.32 ) 

IF(XP. GT.-l.. AND. XP.LE.O.)  XM=10.**(-.185*XP*.375> 
IFIXP.GT.O. .AND.XP.LE. 1 . )  XM=10.**(-.095*XP*.375> 

IF(XP.GT.l.. AND.XP.LE. 2.)  XN»10.**(-.05*XP+.33> 

IF(XP.GT .2. )  XM=10.**(-.03*XP+.29) 

LM( J , I ) =2.*PI*ST*XM/ (RHQGt J )*V ( I )**2. > 

MDR= RHOL ( J ) *P I *LM ( J , I ) **3 . / 6 . 

ADR=PI*LM( J , I )**2. 

MREMO ( J , I ) 3 4 . *P I **2 . *DB* ( RHOG ( J  > / RHOL ( J ) )**. 5*RH0L ( J )*F 
1  *V<I)*L(J)*DT 

MREMA( J , I )=MREMA< J , 1-1 ) -MREMO< J , I ) 

IF(MREMA(J,I).LE.O.)  THEN 
MREMO ( J , I ) =MREMA ( J , I - 1 ) 

MREMA( J , I ) =0. 

L(  J ) =0. 

LL( J ) =1 
END  IF 

MNETT  *MNETT+MREMA ( J , I ) 

LNETT =LNETT +L  ( J ) 

NDR  ( J , I ) =MREMO ( J , I ) /MDR 
SDR  ( J  ,  I )  =NDR  ( J  ,  I  >  *ADR 
NDRTT =NDRTT+NDR ( J , I ) 

SDRTT =SDRTT +SDR ( J , I ) 

PRINT*, ' ' 

PRINT*,’  TIME  OF  MASS  INJECTION*  \T(J> 

PRINT*,’  MASS  IN  DROP  FORM*  ',HREMO(J,I> 

PRINT*,'  DROP  DIAMETER*  ',LM(J,I> 

PRINT* , ’  NUMBER  OF  DROPS*  ' , NDR ( J , I ) 
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O 

■J 


Paje 


C 


c 

c 

c 

c 

c 

c 

c 

c 

c 


PRINT*, ’  SURFACE  AREA  OF  DROPS* 
END  IF 

140  CONTINUE 

MNET ( I ) *MNETT 
LNET  ( I ) *LNETT 
MINNET*MINNET+MIN< I ) 

MATOMI I ) =MINNET-MNET ( I ) 

SDRT  ( I )  *SDRTT 
NDRT  (  I )  *NDRTT 
NDRNET =NDRNET +NDRT ( I ) 

SDRNET = SDRNET  +SDRT ( I ) 

PRINT*, ' '  \ 

PRINT*,’  NUMBER  OF  DROPS  AT  THIS  TIME* 
PRINT*,’  SURFACE  AREA  AT  THIS  TIME* 
PRINT*, ’ ’ 

PRINT*, ’  TOTAL  MASS  ATOMIZED* 

PRINT*,’  TOTAL  NUMBER  OF  DROPS* 
PRINT*,'  TOTAL  SURFACE  AREA* 

PRINT*,' • 

PRINT*, V 
130  CONTINUE 


'  , SDR ( J ,  I ) 


' , NDRT ( I ) 

' , SDRT ( I ) 

’  , MATOMI I) 
’ , NDRNET 
’ , SDRNET 


THIS  PART  OF  THE  PROGRAM  WILL  CALCULATE  THE  BURNING  RATE  AND  THE 
MASS  INVENTORY  AS  IF  THRE  IS  ONLY  BURNING  OCCUR ING. 

500  CONTINUE 
PRINT*," 

C  PRINT*,'  *****•««*«**  BURNING  CALCULATIONS  **************’ 

DO  130  1*1, N 
LBURN ( I ) *  1 

180  BRATEt I ) * 1 . 64*PRESS ( I >**.103 
MINNET*0. 

DO  182  1*1, N 

C  PRINT*, ’ TIME* ’,T( I) 

BAREA< I )*PI*V( I )*DB*DT 
MLEFTBl I,I-1)*MIN(I) 

TMBURN*0. 

TMLEFB*0. 

DO  184  J  *1 , 1 

IF ( LBURN (J ) .NE.O)  THEN 

MBURN ( J , I ) * B AREA ( J ) *BRATE ( I ) *RHOL ( J )*DT 
MLEFTB( J , I )*MLEFTB( J , 1-1 ) -MBURN( J , I ) 

IF(MLEFTB( J,I) .LE.O. )  THEN 
MBURN ( J , I ) *MLEFTB ( J , I - 1 ) 

MLEFTBI J , I ) *0. 

LBURN ( J ) *0 
END  IF 

TMBURN*TMBURN+MBURN ( J , I ) 

TMLEFB=TMLEFB+MLEFTB ( J , I ) 

END  IF 

134  CONTINUE 

MINNET  =MINNET+MIN( I ) 

MINTT ( I ) =MINNET 
TTMLEB ( I ) *  TMLEFB 
TTMBURI I )  *MINNET-TMLEFL' 


% 
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PRINT*  1  1 

PRINT*,’ 'TOTAL  MASS  BURNED  AT  THIS  TIME  '  ,TTMBUR(  I ) 
PRINT*, 'TOTAL  MASS  NOT  BURNED  AT  THIS  TIME' ,TTMLEB( I ) 
182  CONTINUE 

PRINT  OUT  SUMMARY  DATA 

WRITE(20,188> 

WRITE (20, 186) ( T  < I ) ,MATOM(I) ,MNET ( I ) , LNET ( I ) , TTMBUR ( I ) , 
1  TTMLEB(I) ,MINTT(I) ,I=1,N) 

188  FORMAT ( //// ) 

186  FORMAT (2X,F4. 2, 2X,F7 . 4,2X,F7.4,2X,F7.4,2X,F10.7,2X, 

1  F7 . 4 , 2X , F7 . 4 ) 

CLOSE (10) 

CLOSE (20) 

END 
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P»9»  1 


C  ***************************************************************************** 

C 

C  PROGRAM  TO  CALCULATE  THE  AMOUNT  OF  LIQUID  THAT  IS  ATOMIZED  IN  AN  ANNULUS 
C  PROGRAM  NAME  IS  JETC.FTN. THIS  VERSION  THIS  VERSION  USES  ODE  SOLVER  "EPISODE" 
C  WRITTEN  BY  NELSON  MACKEN  31  DEC  86 

C  *****************Hi-*#***************HHt************************^*************** 


J  «******##******###»4HHHI’*»#***»***»***#*»*#**#***#*#*****#»**#'IHt***#*****##«*«« 

c 

c  *******************  MAIN  PROGRAM  SETS  UP  ODE  SOLVER  #**■***********#**#****♦* 

C 

~  ***HI^****#^**1Ht****^****************************#**********tl-**************#**-** 


DEFINITION  OF  VARIABLES 


C 

c 


c 

r 

C 


Y(  1 )  =V( I ) 
Y(2)=AV!I) 

Y (3) =RHOG( I ) 
Y ( 4 ) =RHOL ( I ) 
Y ( 5) =MIN( I ) 

Y ( 6 ) =MATOM 
Y ( 7) =3URF 


VELOCITY  OF  JET 
AREA  OF  ORIFICE 

DENSITY  OF  GAS  IN  COMBUSTION  CHAMBER 
DENSITY  OF  LIQUID  ENTERING  COMBUSTION  CHAMBER 
TOTAL  MASS  THAT  HAS  ENTERED  COMBUSTION  CHAMBER 
TOTAL  MASS  THAT  HAS  BEEN  ATOMIZED 

TOTAL  SURFACE  AREA  OF  DROPLETS  THAT  HAVE  BEEN  ATOMIZED 


DIMENSION  T(60) ,PRESS(60) , Y( 10) ,VP< 60) 
COMMON/ V AR / V I ( 60 ) , AV ( 60 ) , RHOL !  60 ) , RHOG ( 60 ) 
COMMON/ CONST/ I, DT 
COMMON/ OUTPUT /AMLEF , ALEN 


OPEN! 10 ,FILE=' CONST. JETC* , STATUS* ' OLD ' ) 
OPEN (20, FILE = 1  OUTPUT. JETC ' , STATUS= ' NEW ' ) 


N=51 

READ( 10 , 100) (T( I ) , VI ( I ) , VP! I ) , AV! I ) ,RHOL( I ) ,RHOG( I ) , 
1  PRESS! I), 1=1, N) 

WRITE! 20, 100) (T< I ) ,VI ( 1 ) ,VP( I ) ,AV( I ) ,RHOL( I ) , 

1  RHOG! I), PRESS! I) ,1=1, N) 


C  PARAMETERS  FOR  DRIVE 
NI=7 

H0*1 .E-10 
EP3= 1 .E-4 
IERR0R=3 
MF=10 


INDEX* 1 

C  INITIALIZE  VARIABLES 
Y( l )=VI ( 1 ) 

Y(2)=AV( 1 ) 

Y(3)=RH0G( 1 ) 

Y(4) -RHOL! 1 ) 

Y( 5) =0. 

Y(6)=0. 

Y  <  7 )  =0 . 
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f| 

1 
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X0=6. IE-3 
WR I TE ( 20 , 300 ) 

DT=.05E-3 
XOUT =XO+DT 

C  SET  UP  LOOP  TO  SOLVE 
N=51 

DO  2  1=2, N 

CALL  DRIVE(NI,XO, HO, Y, XOUT.EPS, IERROR.MF, INDEX) 

IF( INDEX. NE.O)  STOP 

WR I TE ( 20 , 200 ) T ( I ) , AMLEF , ALEN ,  Y  ( 6 ) ,Y(7) 

XOUT  =XGUT+DT 
2  CONTINUE 
CLOSE < 10) 

CLOSE (20) 

C 

100  FORMAT! 7 (F7.2,2X) ) 

200  FORMAT ( 10X,F4.2,5X,F6.4,5X,F6.4,5X,F7.4,5X,F9.2) 

300  FORMAT ( // ) 

C 

END 

END  OF  MAIN  PROGRAM 

*******#******#********#***#»*»************~IH»**»**ftft»***»****»*«**'lt** 

**  SUBROUTINE  VARC  COMPUTES  PROPERTIES  AT  INTERMEDIATE  TIME  STEPS  ** 

**************#****»*****-**lr  ****#**##*■*■*******#***■***#*•»•*##*■  «•«■*****•** 

SUBROUTINE  VARC ( YD 1 ,  YD2 ,  YD3 ,  YD4 ) 

COMMON/ VAR/ V I ( 60 ) , AV ( 60  > , RHOL ( 60 ) , RHOG ( 60 ) 

COMMON/CONST/ I, DT 

DEFINITIONS  OF  DERIVATIVES 

YD1=  D(VELOCITY)/DT 
YD2=  D<  AREA ) /DT 
YD3=  D(RHOG) /DT 
YD4=  D( RHOL) /DT 

YD1=(VI(I)-VI(I-1))/DT 
YD2= ( AV( I ) -AV( 1-1 ) ) /DT 
YDS  = ( RHOG ( I ) -RHOG ( I - 1 ) ) / DT 
YD4= ( RHOL( I ) -RHOL( I - 1 ) )  'DT 

RETURN 
END 

END  OF  SUBROUTINE  VARC 

***#*******#**#****♦#****#**#***********#*****#*#*##**■*#■**** 

#*  SUBROUTINE  ATOM  COMPUTES  PARAMETERS  FOR  ATOMIZATION  **** 


l 
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C  *******'**-*#****-*****-*#*-»HHHHHf****-**#*#**#*<HH»*-****#********* 

c 

SUBROUT I NE  ATOM ( V , AV , RHOG , RHOL , Y5 , Y6 , YD5 , YD6 , YD7 ) 

C 

COMMON/ OUTPUT/AMLEF , ALEN 
PEAL  MDR , LM , MLEFT , LEN 

DEFINITIONS  OF  VARIABLES 

C 

C  V  =  JET  VELOCITY. CM/SEC 

C  A V=  VENT (ORIFICE)  AREA,CM**2 

C  RHGG=  GAS  DENSITY, GM/CM**3 

C  RHOL=  LIQUID  DENSITY,GM/CM**3 

C  FL IN=  TOTAL  MASS  THAT  HAS  ENTERED  COMBUSTION  CHAMBER, GM 
C  PI=  VALUE  OF  "PI" 

C  ST=  SURFACE  TENSI0N,GM/SEC**2 

C  VISL=  LIQUID  VISCOSITY, GM/CM-SEC 

C  DB=  DIAMETER  OF  BOLT, CM 

C  B , XP , F=  DROPLET  FORMATION  PARAMETERS 

C  LM=  WAVELENGTH  OF  DETACHED  FLU ID, DIAMETER  OF  DROPLET, CM 
C  MDR=  MASS  OF  DROPLET, GM 

C  ADR=  AREA  OF  DROPLET , CM**2 

C  MLEFT , AMLEF  =  MASS  REMAINING  IN  LIQUID  CORE,GM 

C  LEN, ALEN=  LENGTH  OF  LIQUID  CORE, CM 
C  R1,R2=  RADII  OF  FRUSTRUM  FORMED  BY  LIQUID  CORE, CM 
C  SURFA=  OUTER  SURFACE  OF  FRUSTRUM  (LIQUID  CORE) 

C  Y(5) ,Y5=  TOTAL  MASS  THAT  HAS  ENTERED  THE  COMBUSTION  CHAMBER, GM 
Y(6) , Y6=  TOTAL  MASS  THAT  HAS  BEEN  ATOMIZED, GM 
Y<7)=  TOTAL  SURFACE  AREA  OF  ATOMIZED  LIQUID, CM**2 
C  YD5=  MASS  FLOW  INTO  THE  COMBUSTION  CHAMBER , GM/SEC 
r  YD6=  YDOT (6) 

C  YD7=  YDOT ( 7 ) 

C 

C  CONSTANTS 

C 

PI=3. 1416 
ST=71 .6 
VISL= . 08935 
DB*1 .7 

C 

C  CALCULATE  THE  MASS  FLOW  INTO  THE  COMBUSTION  CHAMBER 

C: 

YD5=RH0L*AV*V 

C 

C  VARIABLES  FOR  THE  CALCULATION  OF  THE  MASS  ATOMIZED 

C 

B= ( RHOL / RHOG ) * ( ST / ( V I SL* V ) ) **Z . 

XP=LOG10(B) 

IF(XP.LE.-5. )  F=.01*XP+.09 

IF(XP.GT.-5. .AND.XP.LE. -3. )  F=.02*XP+.14 
IF(XP.GT. -3. .AND.XP.LE. -1. )  F=.0375*XP+. 1925 
IF(XP.GT. -1 . .AND.XP.LE. 1 . )  F=.045*XP+.2 

IF ( XP . GT . 1 )  F= .01 75*XP+ . 2275 

IF (XP.LE. -2. )  XM= 1 0 . ** ( - . 29*XP+ . 22 ) 

IF ( XP.GT . -2. .AND.XP.LE.-l. )  XM=10.#*(-.24*XP+.32) 
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Print  file  "JETC.FTN 


IF(XP.GT.-1 . . AND.XP.LE.O. )  XM=10.**<-. 185*XP+.375> 

IF(XP.GT.O. .AND.XP.LE. 1 . )  XM=10.**(-.095*XP+.375) 

IF(XP.GT. 1 . .AND.XP.LE.2. )  XM=10.**(-.05*XP+.33> 

IF(XP.GT.2.)  XM=10.**(-.03*XP+.29) 

LM=2 .  *P  I  *ST*XM/  ( RH0G*V**2 . ) 

MDR=RH0L*PI*LM**3 . /6 . 

ADR=PI*LM**2. 

CALCULATE  THE  SURFACE  AREA  REMAINING  THAT  IS  NOT  ATOMIZED 
ASSUMING  MASS  LEFT  FORMS  A  "HOLLOW"  FRUSTRUM  WITH  CURRENT 
VALUES  FOR  THE  LIQUID  DENSITY. 

Rl=DB/2.+AV/ (PI*DB) 

R2=DB/2. 

MLEFT=Y5-Y6 

LEN=MLEFT/ ( RHOL*P I * ( 1 . /3 . * ( R 1 **2 . +R 1 *R2+R£**2 . ) - 
1R2**2. ) ) 

SURFA=PI*(R1+R2)*(LEN**2.+<R1-R£)**£.)#*.5 
CALCULATE  MASS  ATOMIZED  PER  UNIT  TIME 


YD6=4 . *PI*RHOL  *V* ( RHOG/RHOL ) ** . 5*F*SURFA 
CALCULATE  SURFACE  AREA  INCREASE  PER  UNIT  TIME 
YD7=YD6*ADR/MDR 

SAVE  MASS  LEFT  AND  JET  LENGTH  FOR  OUTPUT 

AMLEF=MLEFT 

ALEN=LEN 

RETURN 

END 

END  OF  SUBROUTINE  ATOM 

**#*#**#-h*******************  **************************** 


***  SUBROUTINE  DIFFUN  IS  CALLED  BY  ODE  SOLVER  ********* 

***************************1HHt*****************iHH*-**##* 

SUBROUT I NE  D I FFUN ( N , XOUT , Y , YDOT ) 

C. 

DIMENSION  Y(N) ,YDOT (N) 

C:  GATHER  VALUES  OF  VELOCITY,  AREA,  DENS  I TYS 

CALL  VARC ( YDOT ( 1 ) , YDOT ( 2 ) , YDOT ( 3 ) , YDOT ( 4 ) ) 

C  GATHER  THE  MASS  ATOMIZED  AND  THE  SURFACE  AREA  ATOMIZED 
C 

CALL  ATOM(Y( 1 ) ,Y(2) ,Y(3) ,Y(4) ,Y(5) ,Y(6) ,YD0T(5) ,YD0T<6) 
1 , YDOT ( 7 ) ) 


o  o  o  o 


Print  file  ,,JETC.FTN,,  Pi9e  3 


RETURN 
END 

END  OF  SUBROUTINE  DIFFUN 

**************************************************************** 
C.  ******  END  OF  PROGRAM  JETC  -  ODE  SOLVER  TO  BE  BINDED  TO  RUN  **** 

r 

C  **************************************************************** 
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